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ABSTRACT 

We investigate the influence of environment on brightest cluster galaxy (BCG) evolution using a sample 
of 63 clusters at 0.3 < z < 0.9 drawn primarily from the Las Campanas Distant Cluster Survey and 
follow-up V, I, and K' photometry. The luminosity evolution of the entire BCG sample is not adequately 
described by a single evolutionary model. Using the integrated light from the cluster detection as a 
proxy for cluster L x and the suggestion by Burke, Collins, & Mann, we set L x — 2 x 10 44 ergs s -1 
to be the division between high and low luminosity clusters. At high redshift (z > 0.6) BCGs from 
low-La: clusters are fainter, on average, than those from high-L^ clusters and are best modeled as having 
constant luminosity with redshift. The BCGs from high-L^ cluster are best modeled as having a stellar 
population that formed at large redshift (zf orrn > 5) and is passively evolving. However, for the entire 
BCG population, the observed V — I and I—K' colors are well described by a single evolutionary model 
in which the stellar populations have Zf orm > 5 and subsequently passively evolve. We conclude that 
accretion is proportionally more significant for BCGs in lower mass clusters at these redshifts (factor of 
2 to 4 increase in mass since z ~ 1 for the low L x systems; Aragon-Salamanca et al. ) and that the 
accreted matter is in the form of systems with evolved stellar populations. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: formation — galaxies: 
elliptical — surveys 



1. INTRODUCTION 

In standard hierarchical formation scenarios, structure 
formation occurs first in regions of highest overdensity. As 
a result of this accelerated formation, it is generally be- 
lieved, for example, that cluster galaxies are on average 
older than field galaxies. However, studies comparing the 
evolution of galaxies in different environments are diffi- 
cult because identifying equivalent field and cluster pop- 
ulations is an ill-defined problem and because the effects 
of various physical processes (interactions, ram pressure 
stripping, tidally triggered star formation, harassment by 
the global cluster potential, etc.) are poorly understood. 
Selection by morphology, luminosity, color, environment, 
or any other observable property of galaxies has the po- 
tential to bias the sample. 

One galaxy population that can be unambiguously iden- 
tified in a variety of environments, by definition, is the 
dominant galaxy, referred to as the brightest cluster galaxy 
(BCG) in clusters and the brightest group galaxy (BGG) 
in groups. Although the particular galaxy that is the BCG 
or BGG may change with time in any environment, that 
process is part of the evolution one aims to study and 
the sample definition is relatively unambiguous and well- 
suited for comparison to simulations. Over the last 50 
years, many have studied the properties of BCGs from red- 



shifts to ~ 1 (some highlights include Humason, Mayall, 
& Sandage 1956; Sandage, Kristian, & Westphall 1976; 
Hoessel 1980; Schombert 1986; Graham et al. 1996); ini- 
tially with the intent to use BCGs as standard candles, 
later to study galaxy evolution. We revisit the topic of 
BCG evolution in light of recent claims of significant ac- 
cretion (Aragon-Salamanca, Baugh, & Kauffmann 1998; 
hereafter ABK) and environmentally induced differences 
(Burke, Collins, & Mann 2000; hereafter BCM) using a 
large sample of BCGs drawn from the Las Campanas Dis- 
tant Cluster Survey (LCDCS; Gonzalez et al. 2001a). Al- 
though BCG evolution may be quite distinct from galaxy 
evolution in general, its study may shed light on the pro- 
cess of hierarchical structure formation in dense environ- 
ments. 

BCGs are the most massive and luminous galaxies. 
However, the findings that at least some BCGs are not 
drawn from the luminosity function of cluster ellipticals 
(Tremaine & Richstone 1977; Dressier 1978) and that 
BCGs may have different luminosity profiles than other gi- 
ant ellipticals (Oemler 1976; Schombert 1986) suggest that 
they have a different evolutionary history than standard 
cluster ellipticals. Aragon-Salamanca et al. (1993; here- 
after A93) observed BCGs at high redshift (0.5 < z < 0.9) 
in the fC-band and found that while the BCG colors agree 
with the prediction of passive evolution models with high 
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formation redshifts (zf orm J> 2), similar to the results for 
cluster ellipticals (Bower, Lucey, & Ellis 1992), their lu- 
minosities match the prediction of no-evolution models. 
Because stellar populations dim as they age, ABK inter- 
pret this lack of even passive luminosity evolution to imply 
significant, relatively recent mass accretion by the BCGs. 
Parameterizing the necessary mass accretion, ABK esti- 
mate that BCGs have grown by a factor of two to four 
in mass since z ~ 1, a result which they find to be con- 
sistent with the expectation from hierarchical models as 
described by semi-analytic models of galaxy evolution (cf. 
Kauffmann, White & Guiderdoni 1993) 

This rather attractive agreement between observations 
and simulations for the evolution of BCGs over ~ 75% 
of the age of the Universe has recently been questioned. 
The focus of the current discussion is whether all BCGs 
are similar, and in particular whether environment plays a 
critical role in determining the properties of BCGs (Hud- 
son & Ebeling 1997; Collins & Mann 1998, hereafter CM; 
and BCM). BCM study the BCGs of 76 x-ray selected 
clusters in the if-band and find that BCGs from low-La, 
clusters (L x < 2.3 x 10 44 ergs s -1 ) display no luminosity 
evolution. (We will refer to the BCGs from high and low 
h x clusters as high-La and low-La BCGs.) Furthermore, 
they find that the ABK sample primarily consists of low- 
La, BCGs at z > 0.5. Therefore, the BCM observation of 
no luminosity evolution in the BCGs from low-La; clusters 
is entirely consistent with that of ABK. However, they find 
that the luminosities of high-La BCGs have evolved since 
z ~ 1. The stellar populations of their high-La- BCGs are 
consistent with the luminosity predictions of passive evo- 
lution models with a high formation redshift (zf orm J> 5). 
Because high-La, clusters best match the cluster selection 
adopted for comparison to the semi-analytic models, they 
reinvestigate the parameterization of the mass evolution of 
high-L x BCGs and find substantially lower rates of mass 
accretion, at most a factor of ^2 increase in mass since 
z ~ 1. This result suggests that BCGs from high-La clus- 
ters experience a different evolutionary scenario than those 
from low-La clusters. 

Because BCG properties appear to depend upon prop- 
erties of the host cluster, it is imperative that samples are 
drawn from large surveys that span as wide a range of clus- 
ter properties as possible to probe that dependence, and 
that the samples are sufficiently large to average over out- 
liers. Both ABK and BCM have relatively small samples 
at z > 0.5 (~a dozen each) that are biased with respect 
to cluster L x (ABK toward low-La systems and BCM to- 
ward high-La systems). Obviously, enlarging the samples 
of these high redshift clusters has been difficult and any 
potential advance should be exploited. We examine the 
evolutionary history of BCGs using a new sample of high 
redshift clusters drawn primarily from the LCDCS (Gon- 
zalez et al. 2001a). 

The LCDCS is a drift scan survey, conducted by the 
authors at the Las Campanas Observatory in 1995, that 
covers 130 deg 2 of the southern sky and yielded over 1000 
cluster candidates. Rather than selecting overdensities of 
resolved galaxies, our cluster-finding technique relies upon 
the assumption that the total light from a distant cluster 
is dominated by the contribution from unresolved sources. 
In an intrinsically uniform image, the integrated diffuse 



cluster light will manifest itself as a low surface brightness 
enhancement. Not only does this catalog result in a fac- 
tor of 5 increase in the number of cataloged clusters at 
these redshifts, but because our cluster identification cri- 
teria differs from those utilized in previous surveys, this 
catalog provides an independent, well-defined sample with 
which to compare the results from more traditional sur- 
veys. Details of the survey and the catalog of clusters can 
be found in Gonzalez et al. (2001a). 

Because extensive follow-up observations are pro- 
hibitively time-intensive, the BCGs we examine here are a 
small subset of the final catalog. We use 54 clusters from 
the LCDCS and 9 northern clusters from a smaller Palo- 
mar 5m survey (for which the data was originally taken as 
part of the Palomar Transit Grism Survey for high redshift 
quasars; Schneider et al. 1994 and analyzed in Dalcanton 

1995) which we originally utilized to test the low surface 
brightness technique. The combined sample consists of 63 
clusters, 17 of which are spectroscopically confirmed and 
46 of which have photometrically estimated redshifts (Nel- 
son et al. 2001a). We obtain deep follow-up imaging in V, 
I, and K' , although not every BCG is imaged in every 
band. In §2 we describe the data. We include a brief 
description of the survey, cluster candidate selection and 
classification in §2.1. In §2.2 we describe follow-up obser- 
vations and photometry. Next, we give a summary of our 
photometric redshift estimators and BCG selection in §2.3 
and §2.4, respectively. In §2.5 we describe the technique 
used to estimate La;. In §3 we examine the evolution of 
BCGs in both luminosity and color. Finally, we summa- 
rize our results and conclusions in §4. 

2. THE DATA 

Our data originate from a variety of telescopes and in- 
struments. We identify the candidate galaxy clusters using 
drift scans and techniques described briefly below for con- 
text, but in full detail by Gonzalez et al. (2001a). The 
cluster sample and observations used here stem from deep 
optical and infrared follow-up imaging of a small subset of 
the full catalog which was obtained to aid in the classifi- 
cation of candidates and to develop photometric redshift 
indicators. We will briefly review the classification cri- 
teria based upon the deep imaging and our photometric 
redshift estimators. Nelson et al. (2001a) contains a full 
description. 

2.1. Cluster Sample 

We draw our clusters from two optical drift scan sur- 
veys. The first is a ~ 17.5 deg 2 Palomar 5m survey using 
the 4-shooter camera (Gunn et al. 1987) and the F555W 
"wide V" filter (Dalcanton et al. 1995, 1997). Dalcanton 
et al. (1995, 1997) made use of existing scans (originally 
taken as part of the Palomar Transit Grism Survey for 
high redshift quasars; Schneider et al. 1994) to verify the 
low surface brightness (LSB) technique for detecting LSB 
objects to the required central surface brightness. This 
survey identified more than 50 northern hemisphere clus- 
ter candidates. The second is a 130 deg 2 survey of the 
southern sky conducted at the Las Campanas Observa- 
tory lm telescope over 10 nights by the authors using the 
Great Circle Camera (Zaritsky, Shectman, & Brcdthauer 

1996) and a wide filter, denoted as W, that covers wave- 
lengths from 4500 A to 7500 A. This survey was specifically 
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designed to study LSB galaxies and high redshift galaxy 
clusters. We stagger scans in declination by half of the 
field-of-view (24 arcmin) so that each section in the sur- 
vey is observed twice, each time with a different region of 
the detector. The sidereal rate limits exposure times to 
~97 sec per scan at the declination of the survey (—10° to 
— 12°), so that the net exposure time is ~195 sec at any 
location. The pixel scale is 0.7 arcsec/pixel with a typi- 
cal seeing of about 1.25 arcsec. The full region surveyed 
measures ~ 85° x 1.5° and is embedded within the area 
covered by the LCRS (Shectman et al. 1996). The result 
of iterative flatficlding, masking, filtering, and object de- 
tection is a catalog of ~ 1000 cluster candidates (Gonzalez 
et al. 2001a provide a full description of the method and 
catalog) with an estimated 30% contamination rate. 

From the drift scans, we select a set of candidates for 
which to obtain deep follow-up imaging in V, I, and K' . 
These deeper images allow us to classify detections and 
train the cluster-finding algorithm. The cluster candi- 
dates presented in this work were confirmed as clusters 
based upon the following criteria: 1) the appearance of 
the background-subtracted /—band luminosity function, 
in particular a well populated and defined bright end of 
the luminosity function; 2) the presence of a prominent red 
envelope in the V — I CM diagrams, and 3) a well-defined 
concentration in the surface density of galaxies, in partic- 
ular a high degree of spatial clustering of red galaxies. A 
cluster candidate that strongly satisfies any of these three 
criteria is considered a viable cluster. We do not require 
successful candidates to satisfy all three criteria for two 
reasons. First, we do not have the necessary data to eval- 
uate every cluster equivalently. Second, we have confirmed 
clusters that do not strongly satisfy all three criteria (for 
example, high redshift clusters generally have less promi- 
nent red envelopes due to the combination of cosmological 
shifting of light out of the F-band and our modest expo- 
sure times; Nelson et al. 2001a). We acknowledge that 
our classification procedure is subjective, but simulations 
demonstrate that random fields and non-cluster detections 
(LSBs and spurious detections) do not yield well-defined 
luminosity functions nor prominent red envelopes (Nelson 
et al. 2001a). Nevertheless, we will discuss the effect of 
significant contamination on our results in §3.1.1. 

Because the original goal of the follow-up imaging was 
to test our cluster selection algorithms, the cluster sub- 
sample is heterogeneous and thus includes both some of 
our best and some of our worst candidates. In addition, 
since our selection procedure evolved during this testing 
process, some of the confirmed candidates presented here 
are not in the final, statistical cluster catalog presented by 
Gonzalez et al. (2001a). We expect no particular biases 
in BCG properties from this procedure, except those that 
correlate with cluster properties. The potentially influ- 
ential cluster properties include concentration and cluster 
mass, which factor into the cluster selection and which are 
present at some level in all cluster surveys. These biases 
are discussed in detail in §2.4. 

2.2. Optical and Infrared Photometry 

We briefly review the data acquisition and reduction 
process. Readers interested in a more detailed description 
are referred to Nelson et al. (2001a). We image our clus- 
ters in the Johnson V and Cousins /-band because these 



bands straddle the Ca II H & K break for galaxies at red- 
shifts 0.3 < z < 0.9, and K' because it probes the older, 
quiescent stellar populations. Because the majority of our 
photometry was completed before the cluster finding algo- 
rithm was fine-tuned (indeed it was necessary to help train 
the algorithm), we imaged some candidates only in /, to 
provide more images of potential clusters at the expense 
of color information. The data were collected at a variety 
of telescopes from 1996 to 1998. Northern clusters were 
generally imaged at the Palomar 1.5 m telescope and the 
Lick 3 m telescope, while southern clusters were generally 
imaged at the Las Campanas Observatory using the 1 m 
and 2.5 m telescopes. Both the optical and infrared data 
were reduced following standard procedures (see Nelson et 
al. 2001a for details). We calibrated the optical data us- 
ing Landolt's standard fields (Landolt 1983, 1992). The 
K' -band data was calibrated using the HST faint near-IR 
standards (Persson et al. 1998). 

Table 1 contains information regarding the observations 
of each cluster. The first column lists the cluster designa- 
tion. RA and DEC (2000.0) are in the next column, where 
the center of the cluster is taken to coincide with the cen- 
troid of the low surface brightness fluctuation. The last 
six columns contain the exposure time and seeing for each 
of the three photometric bands, V, I, and K' . In sum- 
mary, 97% of our clusters are imaged in /, 54% in V, 35% 
have infrared K' data, and 29% have spectroscopy. Some 
clusters were imaged twice in various photometric bands 
and are listed twice in Table 1 . If one data set is clearly of 
superior quality than the other, we use the superior data. 
Otherwise we average the photometric results. 

The correlation between the luminosity of the BCG mea- 
sured within the fixed radius r m , referred to as the metric 
aperture or L m , and the logarithmic slope of the BCG 
curve of growth, a = [d log(L m )/d log(r)]| r=rm (Hoes- 
sel 1980; Lauer & Postman 1994; Postman & Lauer 1995, 
hereafter PL95; Hudson & Ebcling 1997; CM) is often 
used to place BCG luminosities on a standard-candle scale. 
Such an adjustment typically reduces the scatter in the 
distance-magnitude relation by <~ 0.1 — 0.2 mag. Conse- 
quently, BCG photometry is predominantly measured us- 
ing a fixed aperture of projected physical radius (~10 kpc 
locally, and ^25 kpc at high redshift). However, we choose 
neither to measure our BCG photometry with fixed met- 
ric apertures, nor to correct their luminosities using the 
L m — a relation, for two reasons. First, seeing correc- 
tions that must be applied to fixed aperture magnitudes 
are too uncertain (especially at high redshift, because the 
seeing varied dramatically from run to run, 0.8" — 3" in 
the optical and 0.6" — 2" in the infrared). Second, as we 
describe in detail below, the majority of clusters in our 
sample do not have spectroscopic redshifts. Instead we 
rely upon photometric redshift estimators whose associ- 
ated errors (a p hot ~ 0.06) are comparable to the L m — a 
correction. Applying the L m — a relation to BCG photom- 
etry typically changes Mbcg by ^0.5 mag (see Figure 1 in 
PL95; Figure 10 in CM). For the redshift range in which 
the bulk of our clusters lie, the uncertainty in z p hot corre- 
sponds to AMbcg ~0.2-0.4. Despite the potential increase 
in scatter by not using the L m — a relationship, we argue 
that such precision is lost in comparison to simulations (or 
stellar evolution models) that do not model evolution in 
the structural properties of galaxies. Our sample is use- 
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ful for global tests, but has high photometric scatter that 
prevents case-by-case comparisons. 

Instead of an aperture magnitude, we measure magni- 
tudes using SExtractor's (Bcrtin & Arnouts 1996) "best" 
total magnitudes. For uncrowded objects, this magni- 
tude is the automatic aperture magnitude, which has been 
demonstrated to miss less than 6% of the light from a 
galaxy (Bertin & Arnouts 1996). The corrected isophotal 
magnitude, which misses less than 15% of the galaxy's 
light, but is more robust, is used for crowded objects 
(Bertin & Arnouts 1996). These techniques are not neces- 
sarily correct for galaxies with large extended halos, such 
as those BCGs may have, but the low surface brightness 
halo is missed at all redshifts (see Gonzalez et al. (2001a) 
for a demonstration that ~ 30% of the light from a nearby 
cD is missed by standard techniques). Aside from this 
last uncertainty, the internal uncertainties as calculated 
by SExtractor using the detector read noise and gain, and 
the photon statistics of sky and object, dominate. How- 
ever, we also include the uncertainty from the photometric 
solutions in determining the uncertainty in the final mag- 
nitudes. Galaxy magnitudes are corrected for extinction 
using the dust IR emission maps of Schlegel et al. (1998), 
but are not K-corrected. 

2.3. Spectroscopic and Photometric Redshifts 

To study the evolution of the BCGs, we must know their 
relative distances. We either obtain redshifts directly from 
spectroscopy of a number of potential cluster members, or 
indirectly from photometry of cluster members. Because 
the final cluster catalog contains ^1000 candidates, it is 
not feasible to obtain spectroscopy for all our potential 
clusters, nor even for the 65 clusters discussed here. 

We obtained spectra of 19 of our clusters using the Low- 
Resolution Imaging Spectrograph (Oke et al. 1995) on the 
Keck I telescope in 1995 Dec 20-21, 1997 March 14-15, 
1998 April 4-5, and 1999 March 20-21. Due to weather 
and technical problems ~ 40% of the time was usable. The 
deep imaging necessary to construct multi-object masks 
was generally not available so we used a single long slit 
with a 600 line mm -1 grating. The slit was aligned, using 
the guider image, to include as many individual galaxies as 
possible, generally 2 to 3 with m^ < 22, and to lie across 
the position of the LSB feature detected in the scans. Typ- 
ical spectroscopic exposure times were 30 minutes to one 
hour. Details of the reduction and analysis of the spec- 
troscopy can be found in Zaritsky et al. (1997) and Nelson 
et al. (2001a). 

For clusters without spectroscopy, we estimate the red- 
shifts photometrically using two different techniques. First 
we use the /—band luminosity function measured from the 
/-band follow-up observations to obtain zi um - Second, we 
use the location of the "red envelope" in the cluster color- 
magnitude (CM) diagram measured using the V- and I- 
band follow-up observations to obtain z co ;. For details 
of the procedure see Nelson et al. (2001a). Each tech- 
nique has its strengths and weaknesses. The more practi- 
cal zi um is more efficient than z co i because it requires deep 
follow-up imaging in only one band. However, the cluster 
galaxy catalog must be corrected for interloper contamina- 
tion by statistically subtracting the background galaxies, 
which introduces a significant source of uncertainty. On 
the other hand, z co i does not require background subtrac- 



tion because the red ellipticals in the cluster are typically 
the reddest objects in the field, but the technique is less 
efficient because it requires follow-up observations in two 
separate filters. We empirically calibrate both methods 
using a subset of our clusters with spectroscopic redshifts 
and the 19 clusters taken from the literature (A93; Smail 
et al. 1997; Stanford et al. 1998; see Nelson et al. 2001a 
for details). 

Using clusters with 

Zspec w6 find that z\ ura z S p ec has an 
rms scatter about zero of 0.06, which we adopt as our error 
in zium- Similarly, we calculate z co i — z spec and find the 
rms scatter about zero is 0.07, and adopt this value as our 
error in z co i. These values of the rms underestimate the 
rcdshift uncertainties (because the spectroscopic clusters 
are perhaps richer than the average candidate), but our 
conclusions do not depend on having precise redshifts or 
exact rcdshift uncertainties. All clusters with z spec were 
used for these estimates. For clusters without z spec , we 
rely on the two photometric redshift estimations. For our 
full redshift range, if both zi um and z co i exist, we average 
the two determinations. We disregard any clusters with 
zium > 0.7 that do not also have z co i because zi um becomes 
unreliable at high rcdshift; our observations are generally 
too shallow to contain a sufficient number of galaxies to 
produce a well populated luminosity function. This last 
cut removes nine BCGs from the final analysis. The red- 
shift distribution of our BCGs with spectroscopic redshifts 
(filled histogram) and the full sample (open histogram) is 
shown in Figure 1. 

2.4. Estimating Cluster L x 

As has been recently discussed by BCM and will become 
more evident in §3, environment, as classified on the basis 
of L x , is related to the properties of BCGs. Although x- 
ray observations do not yet exist for our clusters, we found 
that the integrated light of the cluster can be used as a 
proxy for x-ray luminosity (Gonzalez et al. 2001b). Gon- 
zalez et al. obtained drift scan data for 17 known clusters 
with x-ray observations (only one of which lies within the 
original survey area). To directly compare to the results 
from previous high redshift BCG studies (CM; BCM), we 
calculate L x for the EMSS passband (0.3—3.5 keV) assum- 
ing H =50 km s _1 Mpc^ 1 , fl m — 1, and f^A = 0. Non- 
EMSS clusters are converted to the 0.3—3.5 keV passband 
using XSPEC assuming a Mekal model (Mewe, Gronen- 
schild and van den Oord 1985; Mewe, Lemen and van den 
Oord 1986; Kaastra 1992) with a neutral hydrogen col- 
umn density, n# = 5 x 10 20 atoms cm~ 2 , and a metal 
abundance, [Fe/H] = 0.3. For clusters without published 
x-ray temperatures, we assume T=5 keV. Because the ma- 
jority of the non-EMSS clusters were observed by ROSAT 
(0.2—2.5 keV passband), our results are only weakly de- 
pendent upon these assumptions. 

In Figure 2 we plot L x vs £(1 + z) 5 , where £ is the peak 
surface brightness of the detection fluctuation in counts 
s _1 arcsec~ 2 . For L x ;> 2 x 10 44 ergs s _1 the correlation 
with S(l + zf is good. Below this threshold E(l + zf 
systematically underestimates L x and the scatter increases 
significantly. The threshold is remarkably similar to that 
which divides the BCG sample of CM into a high-L^ ho- 
mogeneous population and a low-La, heterogeneous pop- 
ulation, suggesting that the breakdown in our relation is 
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perhaps physically motivated. We adopt L x = 2 x 10 
ergs s _1 as the threshold between a high-Lj. and low-L^ 
sample, which corresponds to S(l + z) 5 ~ 0.100 counts 
s _1 arcsec 2 . 

Using the empirically calibrated L x -S(l + z) 5 relation in 
conjunction with the surface brightness limit of the clus- 
ter catalog, we estimate the limiting x-ray luminosity of 
the entire LCDCS as a function of redshift (Gonzalez et 
al. 2001a). At low redshift the LCDCS is sensitive down 
to the level of poor groups, while at z J> 0.6 we are pri- 
marily detecting high-L^ clusters only. In Figure 3 we 
plot £ as a function of redshift and search for a similar 
trend in the cluster subsample examined in this work. For 
z <J 0.6, we detect clusters with a wide range in L x , while 
above this redshift we become increasingly less sensitive 
to the lower mass clusters. Consequently, any correlations 
between BCG properties and cluster mass will result in 
apparent trends with redshift if we do not correct for this 
selection bias. 

2.5. BCG selection 

Selecting the galaxy that is truly the BCG requires care- 
ful consideration of several factors: passband used for se- 
lection, search radius, and contamination. Even at low 
redshift, where extensive spectroscopy is available, BCG 
selection is problematic. For example, PL95, who stud- 
ied 119 BCGs with z < 0.05 (which included a subset of 
Hocsscl, Gunn, & Thuan's (1980) sample) found that be- 
cause their selection criteria differed from that of Hoessel, 
Gunn, & Thuan (1980), they selected a different galaxy as 
the BCG in seven out of 34 clusters. 

Our first consideration is which passband to use for se- 
lection. Although bright cluster ellipticals have remark- 
ably homogeneous colors (A93; Stanford et al. 1995; Lubin 
1996; Stanford et al. 1998; Kodama et al. 1998; Nelson et 
al. 2001a), even a small spread in colors can result in a dif- 
ferent choice of BCG depending on the bandpass. Of the 
three bandpasscs available to us, K' is the optimal because 
it is most sensitive to the older stellar populations through- 
out the redshift interval of our sample. However, because 
of the small ficld-of-view of IR arrays (corresponding to 
150 — 250ft -1 kpc at the cluster redshifts), K' selection 
will miss BCGs that are significantly displaced from the 
cluster center. Instead, we choose to use the /-band for 
selection because we have large ficld-of-view images in the 
/-band of our entire sample. The /-band lies redward of 
the Ca II H & K break for the majority of clusters in our 
sample (z <; 0.8) and consequently is not overly sensitive 
to recent star formation. 

Our second consideration brings us to the greatest diffi- 
culty in selecting high redshift BCGs: foreground contam- 
ination. Unlike for nearby clusters, extensive spectroscopy 
necessary to confirm cluster membership is unavailable for 
the majority of clusters at high redshift. Furthermore, 
PL95 found that although many BCGs are projected close 
to the cluster center {{tbcg) < 70ft- -1 kpc), a few lie at 
large projected distances {tbcg ~1 Mpc) where contam- 
ination is high. Choosing the optimal search radius, so 
as to include as many BCGs as possible while minimizing 
contamination, becomes a crucial issue. Color informa- 
tion is extremely useful in mitigating contamination be- 
cause field contaminants are lower redshift galaxies that 
are bluer than the high-redshift ellipticals. We have col- 



ors (V — I and/or I—K' ) for 60% of our clusters. To 
determine the optimal search radius, we first identify the 
BCG in clusters for which color information is available. 
We then use the radial distribution of these color-selected 
BCGs to determine the optimal search radius with which 
to maximize our chance of finding the proper BCG for clus- 
ters without color information (i.e. the radius at which we 
minimize the combined impact of foreground contaminants 
and failure to include the BCG). 

For the first step in this sequence, color-aided BCG se- 
lection, we set a large search radius because the selection 
is less susceptible to interlopers. PL95 found that 90% 
of their 119 local BCGs (z < 0.05) lie within a projected 
separation of 350ft _1 kpc. Therefore, for our clusters with 
cither V — I and/or I—K' colors available, we consider 
all galaxies within 350ft -1 kpc of the cluster center, where 
the cluster center is defined to be the centroid of the low 
surface brightness detection on the drift scan images. The 
sole modifier to this criteria is that when we are using 
the K' images, our small field-of-view limits our search 
to the central 150 — 250ft -1 kpc of the cluster. Because 
bright cluster ellipticals occupy a very narrow locus in the 
cluster color-magnitude diagram (cf. SED98), we exclude 
BCG candidates whose colors are 0.4 mag bluer than the 
location of the red envelope in cither V — I or I—K' . The 
location of the red envelope in color space is defined to be 
the position of maximum change in the number of galaxies 
as a function of color (see Nelson et al. 2001a for details). 
Of the remaining galaxies, we choose the galaxy with the 
brightest total magnitude in the /-band as the BCG. 

We examine both the distribution of Tbcg and the 
contamination rates as a function of selection radius for 
color-selected BCGs, r^, to find the optimal search ra- 
dius for those clusters without color information, r™°p~ co1 . 
Figure 4 shows the distribution of Tbcg f° r BCGs se- 
lected by color within 350ft -1 kpc of the cluster cen- 
ter. The open histogram shows the radial distribution 
for all color-selected BCGs, while the filled histogram is 
for color-selected BCGs with z spec . Although there is 
a high concentration (~40%) of BCGs very close to the 
cluster center (tbcg < 50ft -1 kpc), there is also a signif- 
icant number of BCGs (<~40%) at larger projected sepa- 
rations {vbcg > 150ft -1 kpc). Next, we vary the search 
radius and estimate contamination rates, on a per clus- 
ter basis. First, we identify the "true" cluster BCG us- 
ing the color-selection criteria outlined above (including 
r sei = 350ft -1 kpc). Then, we choose 10 random cen- 
ters on the image that lie entirely outside of the cluster 
radius (350ft -1 kpc+r"7 -coi ) and select the brightest I- 
band galaxy that lies within r < r™°™~ co/ without using 
color information. Finally, we define the contamination 
rate for this cluster as the fraction of times that a galaxy 
within the search radius, f™°™~ c °', is brighter than the 
"true" (i.e. the original, color-selected ) BCG. We repeat 
this entire process for every cluster using various search 
radii. Figure 5 shows the distribution of contamination 
rates for four values of r™°p~ co1 . For very small search ra- 
dius, f™ e °" _co ' = 50ft -1 kpc, the contamination rates are 
extremely low, ^80% of the clusters have effectively no 
contamination and there are no clusters with a contami- 
nation rate greater than 20%. By the time 7*™°™ _co ' reaches 
200ft -1 kpc, the contamination rates are quite high, 55% 



6 



Revisiting BCG Evolution with the LCDCS 



of the clusters have contamination rates >30%. We also 
calculate (Am) = m^cG — m cont and find that (Am) <~ 1 
mag in V and /. 

Because there is a significant difference between the 
magnitude of the true BCG and that of a foreground 
contaminant, the optimal search radius for non-color- 
selected BCGs must be small enough to keep contamina- 
tion relatively low, yet large enough to include a signifi- 
cant fraction of BCGs. While the contamination rate for 
50h~ 1 kpc is attractive, such a sample would 
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of BCGs from z ~ 0.3-1 and search for differences in their 
evolutionary histories as a function of environment. 

3.1. BCG Luminosity Evolution 

We search for the influence of environment on evolu- 
tion by examining V, I, and K' Hubble diagrams (Figure 
7) and differentiating between BCGs from high-La, clus- 
ters, low-La, clusters, and clusters for which we do not 
have L x estimates 1 . Similar to CM and BCM, we adopt 
L x = 2 x 10 44 ergs s _1 as the threshold between a high-La, 
and low-La, sample. The errors in the photometry are neg- 
ligible, but the uncertainty in the photometric redshifts is 
not. An error in z p hot of 0.07 corresponds to Amj ~ 0.3 
mag at z ~ 0.6. 

We compare the BCG luminosities to spectral synthe- 
sis models using the GISSEL96 model (Bruzual & Char- 
lot 1993) for an elliptical galaxy experiencing cither no- 
evolution or passive evolution. For all models, we as- 
sume an initial 10 7 yr burst of star formation and a 
Salpeter initial mass functions for masses between 0.1 M 
and 100 Mq. We determine the magnitude zero-points for 
the models by using 64 of PL95's local BCGs whose host 
clusters have x-ray luminosities from ROSAT. Only two 
of these clusters have L x > 2 x 10 44 erg s~ 4 in the EMSS 
passband and consequently this is essentially a low-La, nor- 
malization. For each BCG, PL95 provide photometry in 
a series of apertures of increasing angular size. Because 
we use total magnitudes, rather than fixed metric mag- 
nitudes, we use the largest available aperture from PL95. 
We convert their R c photometry to V, I, and K' using 
the color predictions of the GISSEL96 model (Bruzual & 
Chariot 1993) for a passively evolving elliptical galaxy that 
formed at Zf orm = 10. Finally, we use the locus of the con- 
verted low redshift BCGs to normalize the models. 

First, we investigate whether the predictions from a sin- 
gle evolutionary model fit the observed magnitudes of all 
the BCGs, independent of the La, of their host cluster. Be- 
cause at any given redshift there is considerable observa- 
tional scatter in the magnitudes, we also present a binned 
version of the full data set in redshift intervals of Az =0.1 
in V and / (Figure 7 upper right panel and center right 
panel) and Az —0.2 in K' (Figure 7, lower right panel). 
We adopt cr mean of the binned data as our errors and find 
that they are typically smaller than the size of the points 
and are therefore omitted. In the K' -band, the mean BCG 
luminosities are most consistent with the no-evolution pre- 
diction, but systematically fainter than the model predic- 
tions. To determine whether this offset reflects some prob- 
lem with our data, we compare to the published K-h&nd 
photometry of A93 for 19 other BCGs with z = 0.02-0.92 
(Figure 7, lower right panel). The published photometry 
was measured inside a metric aperture of 50 kpc diameter 
assuming TO = l,Sl\ — and Ho =50 km s _1 Mpc~ 4 
and corrected for galactic extinction using the reddening 
maps of Schlegel et al. (1998), but were not K-corrected. 
Whereas our BCGs show higher scatter than those of A93, 
on average our results are consistent. A small change in 
zero-point would bring the no-evolution model in agree- 
ment with both our data and that from A93. Only two of 

1 In order to estimate L x , we must be able to accurately measure the surface brightness of the cluster detection. However, we cannot do this 
for about one-third of the cluster candidates identified in our earliest analysis, a subset of which are used in this work, because those objects 
are near bright stars. Subsequent catalogs (including the final catalog presented by Gonzalez et al. 2001a) arc drawn from analyses using larger 
masked regions. 



-40% of the BCGs. Instead, we choose to 
= 150/i _1 kpc, which contains a significant 
fraction of BCGs (~60%) yet still has relatively low con- 
tamination rates (70% of the clusters have a contamina- 
tion rate <20%). When we miss the BCG, and instead 
recover the second or third rank galaxy, the difference in 
magnitude is on average 0.5 mag (cf. §3.1). 

Although the majority of BCGs are either elliptical or 
cD galaxies (cf. Graham et al. 1996), this is not necessarily 
universal. Because we do not have morphological informa- 
tion for cluster galaxies, our BCGs are not restricted to be 
of a certain type. However, color selection of BCGs could 
bias our results if some BCGs are bluer than our color 
criteria. An inspection of the cluster color-magnitude dia- 
grams shows that for 50% of the clusters, the red-selected 
BCG is simply the brightest galaxy within 350k" 1 kpc, 
regardless of color. Therefore, color does not play a role 
in the selection of 70% of the BCGs in our sample (includ- 
ing the 40% of our sample for which we do not have color 
information) . For the remaining 50% of the color-selected 
BCGs that are not the brightest galaxy within 350/i _1 kpc, 
the brighter galaxy is bluer by > 1 mag in V — I in 65% of 
the cases. Consequently, our sample is not strongly biased 
toward red BCGs. 

In Figure 6 we present /-band postage stamp images of 
the BCGs used in this study (except for 0944+4732 whose 
image did not reproduce well). The images are 1.5' x 1.5' 
and are centered on the BCG. Table 2 lists the properties 
of the BCGs. Column 1 gives the cluster name. The BCG 
redshift and corrected cluster central surface brightness, 
E(l + z) 5 , are listed in Columns 2 and 3. The next three 
columns give the magnitude of the BCG in V, I, and K' 
. The final two columns contain the color of the BCG 
in V — I and I—K' . In summary, 61% of our sample is 
imaged in V, 96% in /, and 35% in K' . We have color 
information (in either V — I or I—K' ) for 63% of the 
BCGs. Finally, we have L x estimates of 61% of the BCGs 
host clusters. 

3. RESULTS 

According to hierarchical clustering models, the char- 
acteristic scale upon which structure forms increases with 
time. Consequently, low redshift clusters of a given mass 
have experienced a different evolutionary history than 
equivalent mass clusters at high redshift. Galaxy evolu- 
tion studies of cluster and field galaxies at z <J 1 have 
supported this idea by providing evidence that environ- 
ment is linked to the properties of galaxies (e.g. ABK; 
CM; BCM). Here we examine the luminosities and colors 
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our K' BCGs are in high-La clusters and the A93 sample 
also consists principally of low-La clusters. Therefore, we 
conclude that the K' luminosities of BCGs from low-La 
clusters are sufficiently homogeneous to be modeled as an 
apparently non-evolving population (confirming A93 and 
ABK's results, but restricting those results to the low-La 
clusters, in accordance with the results of BCM). 

We have a different situation for BCGs observed in V 
and /. The BCGs have mean luminosities that are consis- 
tent with the no-evolution predictions for z <J 0.6 (Figure 
7, upper right panel and center right panel) but are sys- 
tematically brighter than the no-evolution prediction and 
consistent with the passive evolution of an elliptical galaxy 
that formed at Zf orrn = 5-10 at z> 0.6. For these BCGs, 
a single evolutionary model is insufficient to explain the 
observed magnitudes across the entire range of redshift. 
The upper panels of these diagrams illustrate that, un- 
like for the K' sample, these samples contain significant 
fractions of both high and low-La clusters. Because our 
cluster sample is biased toward more massive clusters at 
z ;> 0.6 (cf. Figure 3), it is likely that the high redshift 
clusters without L x estimates are also high-La systems. If 
this supposition is true, it suggests that BCGs from high- 
La clusters follow a different evolutionary scenario than 
those from low-La clusters. We suggest that the inclu- 
sion of both populations is responsible for the apparently 
erratic Hubble diagrams. 

The differences between the two populations can be ex- 
amined in a variety of ways. First, we examine the de- 
viations from a single evolutionary model in the Hubble 
diagrams more closely. The BCG V and / magnitude 
residuals, Am = tubcg — f^modeU f° r a passive evolution 
model with Zf orm = 5 is shown versus redshift in Figures 
8 and 9 for low-La clusters (left panel) and high-La clus- 
ters (right panel). We omit low-z (z < 0.6) clusters for 
which we could not measure E, but we do include clus- 
ters with z > 0.6 for which we could not measure E (filled 
triangles) because at high redshift we expect to only de- 
tect very massive clusters (cf. Figure 3; Gonzalez et al. 
2001b).' The BCG magnitudes (except for the V-band 
high-L.,, sample) show significant scatter, ~ 0.40 mag. Be- 
cause errors in the photometry are negligible, this scatter 
is most likely due to the uncertainties in the photometric 
rcdshifts (at z ~ 0.6, Az = 0.06 corresponds to Ami ~ 0.3 
mag). The filled boxes denote the mean BCG magnitude 
residuals in V and I with error bars that are a mean . One 
high-La cluster (1024—1239) has very discrepant magni- 
tude residuals (>8cr) in V and / and therefore is omitted 
from the determination of the mean magnitude residuals 
for the high-La sample, under the assumption that it rep- 
resents an incorrect identification of the BCG or a spuri- 
ous detection. In the V-band, (Am 1 ™" 1,1 ) = 0.42 ± 0.12 
and (AmJ, ish_L ') = 0.07 ± 0.10 for low-La and high- 
La clusters respectively, while (AmJ )w_Ll ) = 0.58 ± 0.09 
and (Amf sh - L -) = 0.31 ± 0.11 in the /-band. We con- 
clude that BCGs from high-La clusters are on average 
brighter than those from low-La clusters and more closely 
follow the predictions of a passive evolution model with 

z form — 5 — 10. 

Second, we examine the BCG magnitude residuals from 
the model predictions as a function of S(l + z) 5 . In Fig- 
ure 10, we plot Am = m B ca - m model vs. E(l + z) 5 



in V (left panel) and / (right panel) for passive evolution 
with Zf orm = 5. Filled triangles represent high-z (z > 0.6) 
BCGs without measured values of E. Because our cluster 
sample is biased toward very massive systems at z J> 0.6, 
we consider these BCGs to be from high-La clusters and 
assign them lower limiting values of E(l + z) 5 = 0.100. 
Amy correlates inversely with E(l + z) 5 at the 90% con- 
fidence level according to the Spearman rank test with or 
without the upper limit points included. In /, the correla- 
tion is less significant. Including the upper limits, the cor- 
relation is significant with 90% confidence, but excluding 
them the correlation is significant at only the 85% level. 
These correlations illustrate that as E(l + z) 5 increases, 
the residuals from the passive evolution model, tend to 
decrease. 

We summarize the observational situation as follows: 1) 
the BCG sample includes clusters of both high and low 
La, where the dividing line in L x is set, both from previ- 
ous studies and the behavior we observe between L x and 
E, at La = 2 x 10 44 ergs sec^ 1 ; 2) K' observations of BCGs 
(both our data and that of A93) are consistent with no- 
evolution models, but are based on BCGs that are almost 
exclusively in low-La clusters; 3) the V and I BCG Hubble 
diagrams, which include BCGs from both high and low-La 
clusters, cannot be fit with a single class of evolutionary 
model; and 4) after dividing the sample according to host 
cluster La we find that the BCGs from low-La clusters 
fit the no-evolution model in both V and I (in agreement 
with the results from our K' observations), and that the 
BCGs in the high-La clusters are better fit by a passive 
evolution model, particularly in V where recent star for- 
mation should be more detectable. To address whether 
these results are driven by physical effects, wc next focus 
on potential problems with our data and approach. 

3.1.1. Potential Pitfalls 

We observe that the luminosities of high redshift BCGs 
(z ;> 0.6) are, on average, brighter than the no-evolution 
extrapolation of the luminosities of their low redshift coun- 
terparts. Because our sample is biased toward more mas- 
sive systems at high redshift, we interpret this as evidence 
that environment may influence the evolution of BCGs. 
However, there are other factors that could explain the 
observed "brightening" of BCG luminosities at z ;> 0.6. 
We discuss the possibilities below. 

False Clusters: The fraction of false detections of clus- 
ters in the original survey increases as we consider larger 
rcdshifts (cf. Gonzalez et al. 2001a) and one might suspect 
that this effect could lead to various misleading trends in 
the "apparent" BCG population. Because of our follow-up 
imaging, we do not consider contamination to be as serious 
for this cluster subsample as for the survey itself. First, 
our redshift estimators (the luminosity function and red 
envelope method) fail to converge on random (i.e. spuri- 
ous) centers and on non-cluster (low surface brightness 
galaxies) detections in our follow-up images (Nelson et 
al. 2001a). Second, the luminosities of randomly selected 
brightest field galaxies do not follow the predictions of ei- 
ther passive evolution or no-evolution models for BCGs. 
To illustrate this point, we use the simulations outlined in 
§2.5 and select the brightest field galaxy in an area compa- 
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rable to that of each cluster. Assuming the randomly se- 
lected galaxies lie at the cluster redshift, the dispersion in 
their magnitudes is significantly larger, 07 = 1.2 mag, than 
that of our cluster BCGs (Figure 11, upper panel). Their 
magnitudes ((m/) = 19.2 mag) are typically fainter than 
our BCGs, and do not correlate with the cluster redshift 
(Figure 11, lower panel). Indeed, Figure 11 demonstrates 
that the typical luminosity of a bright field galaxy is fainter 
than or consistent with the predictions of the no-evolution 
models for z <J 0.8. Figures 8 and 9 highlight such a poten- 
tial false cluster detection. 1024—1239, a high-L^ BCG at 
z = 0.57, has magnitude residuals in V and / that are >8cr 
fainter than the mean magnitude residuals of other high- 
ly BCGs. From our simulations, 35% of randomly selected 
brightest field galaxies at 0.50 < z < 0.60 have magnitudes 
that are > 8a fainter than the mean magnitude of the 
BCGs in the same redshift range. Consequently, although 
our sample may include a few false cluster detections, an 
increase in the fraction of false clusters with z in unlikely 
to cause significant brightening of our observed BCGs. 

Biased Measurement of E: If S is dominated by 
the BCG halo, then there will be a correlation between 
brighter BCGs and E. Gonzalez et al. (2000) performed 
a detailed analysis of the distribution of luminous matter 
in the nearby galaxy cluster Abell 1651 and found that 
the BCG contributed 30% of the total cluster light within 
r = 500/i _1 kpc of the cluster center. Thus, in this in- 
stance E could be contaminated by as much as 30% by 
the BCG. If the BCG halo contributes excess light to the 
measurement of E then we expect E to decline with in- 
creasing rsca- Although a slight trend appears in Figure 
12, we find no statistically significant correlation between 
E(l + z) 5 and Tbcg- Therefore, although the question of 
E contamination by the BCG halo remains open, we do 
not find evidence in our sample of biased measurements of 
E. 

Poor Photo-z's: A systematic overestimation of the 
cluster photometric rcdshifts at z ^ 0.6 would cause the 
inferred luminosities of the BCGs to be brighter than their 
actual luminosities. We believe this is not a problem for 
two reasons. First, contamination systematically causes an 
underestimation of the redshift because foreground galax- 
ies tend to be brighter than cluster galaxies (see Gonza- 
lez et al. 2001a for an illustration of the magnitude of 
the effect for photo-z's determined from the BCGs). Sec- 
ond, we find no difference in the behavior of the magni- 
tude residuals for BCGs in clusters with photometric and 
spectroscopic redshifts. In Figure 13 we plot the mag- 
nitude residuals in V (left panel), I (center panel), and 
K' (right panel) for a passively evolving elliptical galaxy 
with Zf orm = 5, differentiating between BCGs with z spec 
(filled circles) and those with z p hot (open circles). Re- 
call that we disregard any clusters with zi um > 0.7 that 
do not also have z co i, because zi um becomes unreliable at 
high redshift. There is no statistically significant differ- 
ence between the luminosities of BCGs in clusters with 
spectroscopically and photometrically estimated redshifts. 
We conclude that that errors in our photometric redshift 
estimations are not causing the observed brightening of 
BCGs at z £ 0.6. 



Contamination: Another concern is that the observed 
brightening of BCG luminosities in V and I at z ^> 0.6 is 
due to increased contamination in the BCG sample with 
redshift. Figure 7 shows that for z <J 0.6 an interloping 
field galaxy would typically need mj <J 19 to be brighter 
than the expected BCG. We find that <5% of all the field 
galaxies in our images are brighter than the expected I- 
band magnitude of the BCG at z <^ 0.6. Consequently, we 
expect insignificant contamination at low redshift. How- 
ever, for z ;> 0.6 the number of potential contaminating 
field galaxies approximately doubles (13% of our compos- 
ite field galaxy population have mj < 20 mag). 

To address this issue in a consistent manner, we use 
the simulations outlined in §2.5 to evaluate the contam- 
ination rate of both non-color-selected and color-selected 
BCGs as a function of redshift. We are chiefly interested 
in whether an increasing contamination rate with redshift 
is responsible for the observed brightening of BCGs rel- 
ative to the no-evolution model prediction at z ^ 0.6. 
Therefore, we compare the magnitudes of potential con- 
taminants to the no-evolution model predictions rather 
than directly to our data. We find that the contamination 
rate for non-color-selected BCGs at low redshift (z < 0.5) 
is quite low, ~10%. However as redshift increases and 
we begin to sample fainter into the luminosity function 
of the interloping galaxies, the contamination rate triples, 
reaching ~30% at z > 0.5. On the other hand, the con- 
tamination rate for color-selected BCGs is quite low and 
remains relatively constant with redshift - 5% contamina- 
tion rate at z < 0.5 and 10% at z > 0.5. This constancy 
of contamination rate arises because bright interlopers are 
typically blue foreground galaxies which are more easily 
differentiated as the cluster E/S0 sequence reddens with 
redshift. 

Are these estimated contamination rates at z > 0.6 
high enough to explain our observed brightening? Al- 
though contamination increases to a modest rate of 30% at 
z > 0.5 for non-color-selected BCGs, only ~20% of BCGs 
at z > 0.6 are selected in this manner. The remaining 80% 
of BCGs are selected using color information which has a 
rather low contamination rate of 10%. From our simula- 
tions, we find that at z > 0.6, foreground contaminants 
are on average 0.63 mag brighter than non-color-selected 
BCGs. For color-selected BCGs, contaminants are 0.41 
mag brighter than the mean BCG magnitude. Assuming 
that 30% of non-color-selected BCGs are actually interlop- 
ers that are 0.63 mag brighter than the BCG, while 10% 
of color-selected BCGs are foreground contaminants that 
are 0.41 mag brighter than the BCG, we find that con- 
tamination causes the mean BCG magnitude to be ~0.1 
mag brighter than the no-evolution model prediction at 
z = 0.6 — 0.8. However, we observe an actual brightening 
of ^0.5 mag for this redshift range in our data. Therefore, 
while contamination may be responsible for some of the 
observed brightening, it is insufficient to explain the full 
extent of brightening of BCG magnitudes relative to the 
no-evolution model predictions at z > 0.6. 

A final way to test whether contamination is responsi- 
ble for our observed brightening is to restrict the sample to 
color-selected BCGs because they have a low contamina- 
tion rate (which remains relatively constant with redshift) . 
The V Hubble diagram is already entirely color-selected 
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because all of our BCGs with U-band data also have I- 
band data. The / Hubble diagram, however, does have 
a significant contribution from non-color-selected BCGs. 
We reproduce the I Hubble diagram (Figure 14, upper 
panel) distinguishing between color-selected vs. non-color- 
selected BCGs. In the lower panel of Figure 14 we bin the 
color-selected galaxies only and again find that the aver- 
age BCG luminosity is systematically brighter at z ;> 0.6. 
Therefore, we conclude that our observed brightening of 
BCGs at z ;> 0.6 is not due to increased contamination 
with redshift. 

Missed BCGs: BCGs have been found to lie as far as 
~1 Mpc from the cluster center (PL95). For those cases 
in which the BCG lies outside our search radius, we will 
miss it and instead choose the second- (or lower-) ranked 
cluster galaxy. Locally, ^90% of BCGs are within 350/i _1 
kpc. Adopting this value for higher redshift clusters, we 
expect that only ~10% of our color-selected BCGs are, in 
reality, the second- (or lower-) ranked galaxy. However, 
the radial distribution of our color-selected BCGs (Fig- 
ure 4) suggests that we may be missing as many as 40% 
of the BCGs for clusters that do not have color informa- 
tion (r™"- co! = 150ft- 1 kpc). How docs this affect our 
results? At high redshift (z ;> 0.6) our BCGs are predom- 
inantly (80%) color-selected. Therefore, we estimate that 
we select the true BCG in at least 75% of our high redshift 
clusters. At low redshift, the majority (~65%) of clusters 
are also selected using color information. For the modest 
number of clusters that are non-color-selected, we may be 
missing as many as 40% of the BCGs. However, we find 
that the mean magnitude difference, Am=m 2 „r bibcg, 
is ^0.5 mag for both color-selected and non-color-selected 
second-ranked galaxies regardless of redshift. Because we 
find similar values of Am for the color-selected and non- 
color-selected samples and no trend of Am with redshift, 
we conclude that if we identified the second-ranked (or 
lower) galaxy rather than the BCG in some systems this 
error did not lead to a detectable effect. 

We conclude that systematic effects are unlikely to 
account for the observed "brightening" of high redshift 
(z ;> 0.6) BCGs in a sample with mixed cluster selection 
that favors the selection of high-L^ clusters at higher red- 
shifts, but warn that such problems need to be kept in 
mind in all such studies (even spectroscopically-selected 
BCG samples are susceptible to some of these pitfalls). 

3.2. BCG Color Evolution 

Evolutionary constraints based on luminosities measure 
only one aspect of the models. For example, a galaxy can 
brighten in any particular passband by either having in- 
creased star formation or adding, via accretion, more stars 
to the initial population. Colors enable us to distinguish 
between these two possibilities and so we now turn our 
attention to the colors of BCGs. Are the colors consistent 
with a stellar population that formed early (zf orm > 5) 
and experiences passive evolution (as indicated by the lu- 
minosities of BCGs in high-L^ environments)? Or are they 
consistent with no-evolution and accretion (as indicated by 
the luminosities of the BCGs in low-La; environments)? 

In the left panels of Figures 15 and 16 we compare the 
V — I and I—K' colors of BCGs to the spectral synthe- 



sis predictions of Bruzual & Chariot (1993) for an ellip- 
tical galaxy experiencing cither passive or no-evolution. 
There is considerable observational scatter in the colors of 
these galaxies, so we bin the data in redshift intervals of 
Az = 0.1 for V — I and Az = 0.15 for I—K' (right panels). 
The behavior in both colors corresponds best with passive 
evolution models of an elliptical galaxy that formed at 
Zform <, 10. Unlike the case of luminosity evolution, one 
evolutionary model is sufficient to explain the observed 
colors of all BCGs, regardless of their environment. It is 
reassuring that the colors do not favor the non-physical 
model of no-evolution for the BCGs in low-La, clusters. 
We conclude that environment does not play a role in the 
BCG colors by observing that there is no correlation be- 
tween V - I vs. S(l + zf (Figure 17). 

4. DISCUSSION 

Combining the constraints from luminosities and colors, 
the results suggest environment plays a role in the mass 
evolutionary history of BCGs. The colors of BCGs from 
both high-La; and low-La; clusters over the redshift range 
z = 0.3-1 imply that they are predominantly comprised 
of old stellar systems (zf orm ^5-10) that have undergone 
little, if any, recent star formation. However, tracing the 
luminosities of high-La; and low-La; BCGs as a function 
of time reveals subtle, but important, differences in their 
evolutionary histories. BCGs in x-ray luminous clusters, 
L x > 2 x 10 44 ergs s _1 , follow passive evolution model pre- 
dictions typical of an old quiescent stellar population. In 
addition, these BCGs have luminosities that are system- 
atically brighter than their low-La; counterparts at high 
rcdshifts. BCGs in low-La; clusters do not exhibit the char- 
acteristic dimming associated with the aging of its stellar 
populations. Instead, their luminosities are fainter on av- 
erage at high redshift relative to a model of a passively 
evolving population and remain relatively constant with 
time. 

ABK interpreted the lack of passive luminosity evolution 
of low-La, BCGs as a sign of mass accretion - BCGs ac- 
crete mass in stars to counter the dimming of their under- 
lying stellar populations, thereby maintaining a roughly 
constant integrated luminosity. Using if-band photom- 
etry of 25 BCGs at < z < 1, ABK parameterized 
the luminosity evolution by Lk{z) = Lk(0) x (1 + z) 7 . 
Assuming H — 50 km s 1 Mpc -1 , the least-squares fit 
to the data yielded 7 = —1.2 ± 0.3 (zf orm = 2) and 
7 = -0.7 ± 0.3 (z form = 5) for VL m = 0.0, Vl K = 0.0, while 
7 = -2.2±0.3 {z form = 2) and 7 = -1.7±0.3 (z form = 5) 
for TO = 1, fl\ = 0.0. Provided that -ftT-band light traces 
mass, they conclude that BCGs increased their mass by a 
factor of 2 to 4 since z ~ 1. In such a model, the galaxy 
colors would evolve passively if the objects accreted had a 
similar star formation history as the BCG. Because ABK 
did not have x-ray observations of their clusters, they as- 
sumed that their inferred mass accretion rate was typical 
of all BCGs. However, BCM re-examined ABK's sample 
with x-ray data and found that the majority of their clus- 
ters have La; < 2.3 x 10 44 ergs s _1 at z > 0.5. To assess 
the effect of cluster mass on inferred BCG mass accretions 
rates, BCM studied an x-ray selected sample of 76 BCGs 
at 0.05 < z < 0.83 and parametrized the mass accretion of 
their high-La; BCGs only. Using a parameterization of the 
same form as ABK, they found that 7 = -0.93 ± 0.23 
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{zform = 2) while 7 = —0.38 ± 0.22 (zf orm — 5) for 
n m = 1.0, A = 0.0 and H = 50 km s^ 1 Mpc" 1 . Con- 
sequently, they conclude that the mass of high-L x BCGs 
has increased by a factor of 2, at most. 

Because our data do not span a large enough baseline in 
rcdshift, we cannot parameterize the BCG mass accretion 
rate as a function of rcdshift for this sample. However, 
in Figure 18, we compare our measured /-band luminos- 
ity evolution to the parametrized mass accretion rates of 
ABK and BCM. The filled boxes denote the mean mag- 
nitude residuals of our low-L^ (left panels) and high-L^ 
(right panels) BCGs for a passive evolution model with 
Zform — 5 (solid line) assuming Vt m = 0.0, il\ — 0.0 (up- 
per panels) and Q m = 1,Oa = 0.0 (lower panels) and 
Ho = 50 km s" 1 Mpc -1 . Although the mass accretion 
rates predicted by ABK (dotted line) were derived using 
predominantly low-L x BCGs, we reproduce the relation in 
all four panels for reference. The dashed line is the evo- 
lution parameterization of BCM for high-L^ BCGs with 
Q m — 1 only (lower right panel). Figure 18 (left panels) 
demonstrates that our low-L^ BCGs are in excellent agree- 
ment with the mass accretion rates of ABK, implying that 
BCGs from low x-ray luminosity clusters accrete a factor 
of ~ 2 if Q m = 0.0 or ~ 4 if Q m = l in mass since z ~ 1. 
Turning our attention to the high-L x sample (Figure 18, 
right panels), we find that the BCGs from high-L^ clusters 
show less mass accretion than those from low-L^ clusters, 
which is consistent with the findings of BCM. However, 
for the only case with which we can directly compare to 
the results of BCM (Cl m = 1, Q\ = 0), our high-L^ sample 
suggests a somewhat higher accretion rate than BCM (a 
factor of ~ 2 since z ~ 1 as opposed to a factor of ~ 1.3). 
In sum, our data confirm that low-L x BCGs may have ac- 
creted a factor of 2- 4 in mass since z ~ 1, while high-L^ 
BCGs are limited to ^ 2 increase in mass since z <~ 1. 

Although some low-L^, BCGs may be accreting at least 
a factor of 2 in mass since z ~ 1, we find little evidence 
for any recent star formation in the colors of BCGs. Be- 
cause the accretion is probably dominated by the infall 
due to dynamical friction of massive early-type galaxies 
near the BCG, large bursts of star formation are presum- 
ably unlikely. However, even a modest amount of gas in 
the progenitors could induce some star formation. If as 
little as 10% of the merger remnant's mass is converted 
into stars, the galaxy's V — I and I—K' colors get bluer 
by as much as 2.5 magnitudes (Nelson et al. 2001a). Al- 
though the duration of the peak bluing is short, <;130 Myr, 
the galaxy is 0.1 mag bluer than the red envelope for ~1 
Gyr, and we should see that reflected in the BCG colors. 
Furthermore, our BCG color selection criteria does not 
significantly hamper our ability to detect a galaxy under- 
going a burst of star formation. A starbursting galaxy is > 
0.4 magnitudes bluer than the location of the red envelope 
for only ^260 Myr, which translates to ~35% of the dura- 
tion of the bluing period (i.e. 0.1 mag bluer than the red 
envelope). From our data, we can only make statistically 
general claims. Even though some BCGs may appear bluer 
than the models, we conclude on the basis of the average 
agreement with the colors from passive-evolution models 
that the dominant form of accretion by BCGs since z ~ 1 
consists of gas-poor, old stellar systems. 

This discussion of BCG evolution suggests that, to some 
extent, observational efforts at studying cluster galaxy 



evolution may be self-defeating. Studies attempting to 
go to higher redshifts (which naturally identify only the 
richest clusters) may face difficulties because although one 
is selecting systems at an earlier time, these systems may 
evolve faster — thereby, countering the expected difference 
between nearby and distant clusters (Kauffmann 1995a, 
1995b). Despite the observations of increased merging 
in some high redshift clusters (van Dokkum et al. 1999), 
galaxy evolution in the observationally accessible regime 
out to z ~ 1 may be most dramatic in poor clusters and 
rich groups (more common in the local samples). Even 
at low redshift, some of these environments show evidence 
for hierarchical evolution (Zabludoff & Mulchaey 1998), 
and the study of stellar populations in the giant ellipticals 
that result in the dynamical collapse of such systems (cf. 
NGC 1132; Mulchaey & Zabludoff 1999) may provide a 
nearby example of the processes that formed the massive 
ellipticals in the richest systems at z ~ 1. 

5. SUMMARY 

We investigate the influence of environment on BCG 
evolution using a sample of 63 clusters at 0.3 < z < 0.9 
drawn from the Las Campanas Distant Cluster Survey. 
This sample represents a factor of 2 increase in the number 
of BCGs studied at z > 0.5. Using the peak surface bright- 
ness of the cluster detection fluctuation, S, as a proxy for 
cluster L x (Gonzalez et al. 2001b), we define L x = 2 x 10 44 
ergs s _1 as the threshold between a high-L^ and low-L^ 
sample. We use this division to trace the luminosity and 
color evolution of BCGs in V, I, and K 1 as a function of 
environment. Our primary results are: 

1) We find that the luminosity evolution of BCGs from 
a mixed sample of high-L^ and low-Lj, clusters cannot be 
adequately explained by a single evolutionary model. 

2) We find that the V, I, and K' luminosities of BCGs 
from low-Lj, clusters are fainter, on average, than the lu- 
minosities of high-Lj, BCGS. The low-L^ BCG luminosi- 
ties are most consistent with no-evolution model predic- 
tions. Our confirmation of the observational results of 
previous work suggests, in accordance with previous es- 
timates (ABK), that these BCGs may be in the process of 
gaining a factor of 2 to 4 in mass via accretion since z ~ 1. 

3) Tracing the luminosities of high-L^ BCGs, we find they 
are better described by passive evolution models with a 
high redshift of formation, Zf orm > 5. Again our agree- 
ment with previous studies (BCM) indicates these galaxies 
are accreting more slowly, a factor of 2 in mass, at most, 
since z ~ 1. 

4) Combining the low-L x and high-L^ results, we conclude 
that we are seeing evidence for the mass-dependent evolu- 
tion rates predicted in hierarchical models. The structures 
in the more massive clusters are closer to reaching their fi- 
nal state, and so are proportionally accreting more slowly, 
while the lower mass clusters are still in a more vigorous 
state of dynamical evolution. 

5) From the colors of BCGs (V - I and I-K' ), we de- 
termine that a single evolutionary model is sufficient to 
describe the observed color evolution of the stellar popu- 
lations of all BCGs. The colors of the combined sample 
of high-L^ and low-Lj, BCGs are most consistent with the 
predictions of passive evolution models with Zf orm > 5. 
We find no evidence in the BCG colors for induced star 
formation and conclude that the dominant form of accre- 
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tion by BCGs since z ~ 1 consists of gas-poor, old 
lar systems (stellar populations similar to those already 
present in the BCG). 

Because the evolutionary history of BCGs appears to be 
linked to properties of the host cluster, it is significant that 
we confirm previously observed trends with a large sample 
of BCGs that have completely different selection criteria. 
Unfortunately, our photometric data only allow us to spec- 
ulate on the merging history of BCGs. Direct evidence of 
accretion events can be obtained only with spectroscopy 
or high resolution imaging. Spectroscopy will yield de- 
tailed information about the current star formation in the 
galaxy and its kinematics but is available for only a limited 
number of galaxies at these redshifts (cf. van Dokkum & 
Franx 1996, Kelson et al. 2000, Kelson et al. 2001). High 
resolution imaging will complement the spectroscopy be- 
cause signatures of merging events include an increase in 
the effective radius and the slope of the surface brightness 
profile of the BCGs (Nelson et al. 2001b, van Dokkum & 
Franx 2001). Tracing the structural parameters with red- 
shift will give an indication of the accretion rate of BCGs. 
The combination of a large cluster sample, such as the one 
presented here, and spectroscopy or high resolution imag- 
ing has the potential to directly address the question of 
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Fig. 1. — The redshift distribution of our sample of BCGs with spectroscopic rcdshifts (filled histogram) and the full sample (open histogram). 
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E*(l + z) 4 {cnts s 1 arcsec 2 ) 

Fig. 2. — Comparison of L x with S(l + z) 5 for 17 known x-ray clusters for which we have drift scan data. h x is the cluster x-ray luminosity 
in the EMSS band (0.3-3.5 keV) and £ is the peak surface brightness of the detection fluctuation. The line is the least squares fit to the data. 
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Fig. 3. — Comparison of S(l + z) 8 with z for clusters used in this work. Because S(l + z) 5 correlates with L x , our sample is biased toward 
more massive systems at high redshift. 
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Fig. 4. — The distribution of the radial separation of the BCG from the cluster center, rsCG, f° r all color-selected BCGs (open histogram) 
and for those with z sp ec only (filled histogram). 



16 



Revisiting BCG Evolution with the LCDCS 



30 r 



-i 1 r 



r~4 1 i 



J3 



10 - 




30 



!■ . n ! 



<50 h*" 1 kpc 



I 1 I I I I I I I I 



h 1 kpc 



10 - 







1~\ 



r^, nM C&J <100 h" 3 kpc 



I I | I I I I 



_GL 



i_CliL,_L 



0.5 

contamination rate 



10 0.5 1 

contamination rate 



Fig. 5. — The distribution of contamination rates for four searching radii: r™°™ co1 = 50h 1 kpc (upper left panel), r™°™ co1 = 100/j 1 
kpc (upper right panel), r™°™~ co1 = \50h~ 1 kpc (lower left panel), and r™°™~ co1 = 200/i~ 1 kpc (lower right panel). The contamination rate 
is defined as the fraction of times that a galaxy within the search radius is brighter than the color-selected BCG (see text for details). 
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Fig. 6. — Postage stamp images of the BCGs used in this work. All images are in the /-band except 1100+4620 which is imaged in 
K' . The orientation of the images is such that north points up and east to the right. The images are of (starting at the top, from 
left to right): 0915+4738, 0936+4620, 1002-1245, 1002-1247, 1005-1209, 1007-1208, 1012-1243, 1012-1245, 1014-1143, 1015-1132, 1017- 
1128, 1018-1211. (Due to size restrictions the astro-ph version only shows one of the set of five panels. The others are available at 
ngala.as.arizona.edu:/dennis/clusters.html or the published version) 
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Fig. 7. — The V-band (upper panels), /-band (center panels), and if '-band (lower panels) Hubble diagrams for BCGs from high-L^ clusters 
(left panels, filled circles), low-L^ clusters (left panels, open circles), and clusters for which ~L X cannot be estimated (left panels, crosses). 
The threshold between high-L^ and low-L^ is 2 X 10 44 ergs s — 1 . The errors in the photometry are negligible, but an error in z p f lot of 0.07 
corresponds to Am; ~ 0.3 mag at z ~ 0.6. The BCG luminosities are compared to the spectral synthesis model predictions of Bruzual & 
Chariot (1993) for an elliptical galaxy experiencing passive evolution and no-evolution. The low redshift points are BCGs from PL95 for 
which the host cluster has x-ray data and are used to normalize the models. In the right panels, our data are binned in Az =0.1 for the V- 
and /-bands (while Az = 0.01 for PL95). For the K' -band we bin the data in Az = 0.2 and compare to the if -band BCG sample of A93 
(open squares), which is not binned. The error bars, which are o mean of the binned data, are smaller than the points and are omitted. 
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Fig. 8. — BCG V magnitude residuals for low-L^ clusters (left panel) and high-L^ clusters (right panel). Because we detect only very 
massive clusters at high rcdshift, clusters without measured values of £ at z > 0.6 are included in the high-L^ sample (filled triangles). 
Residuals are defined as Am = mgcG ~ m modeli f° r a passive evolution model with Zf orm = 5. The filled boxes are the mean magnitude 
residuals and have error bars that arc a mean . The lines are the same as those in Figure 7. 
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Fig. 9. — Same as for Figure 8 but in the /-band 
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Fig. 10. — Comparison of BCG magnitude residuals, Am = rriBCG ~ m modeh m V (left panel) and / (right panel) with S(l + z)° for 
passive evolution with zj orm = 5. Filled triangles are BCGs at z > 0.6 without measured values of E. Because we are only sensitive to 
massive clusters at high redshift, we assign these BCGs lower limiting values of E(l + z) 5 = 0.100. 




Fig. 11. — / Hubble diagram for the brightest field galaxies selected in random areas comparable to that of our clusters (upper panel). 
the lower panel, the data are binned in Az = 0.1. The error bars, which are a mean of the binned data, are smaller than the points and 
omitted. 
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Fig. 12. — Comparison of S(l + z) 5 with tbcg- The lack of a statistically significant correlation between E(l + z) 5 and tbcg 
that the measurements of £ are not biased by contamination from the BCG halo. 
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Fig. 13. — BCG magnitude residuals in V (left panel), I (center panel), and K' (right panel) for clusters with z spec (filled circles) and 
z phot (open circles). Residuals arc defined as Am = mgcG ~ m model, for a passive evolution model with Zf orm = 5. The lines are the same 
as those in Figure 7. 
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Fig. 14. — The /-band Hubble diagram from Figure 7, differentiating between color-selected and non-color-selected BCGs. Symbols are the 
same as for Figure 7. 
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Fig. 15. — BCG colors in V — I compared to the spectral synthesis model predictions of Bruzual & Chariot (1993) for an elliptical galaxy 
experiencing passive evolution and no-evolution (left panel). In the right panel, the data are binned in Az = 0.1 and have error bars which 
are cr mean of the binned data. 
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Fig. 16. — BCG colors in I—K' compared to the spectral synthesis model predictions of Bruzual & Chariot (1993) for an elliptical galaxy 
experiencing passive evolution and no-evolution (left panel). In the right panel, the data are binned in Az = 0.15 and have error bars which 
are ffmeon of the binned data. 
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Fig. 17. — Comparison of the colors of BCGs in V — I with S(l + z) 5 . There is no statistically significant correlation between the BCG 
colors and E(l + z) 5 . 



Nelson et al 29 



2 



2 

II 1 

ijH 

N 

l-n 

E 

<r 

o 



II 1 1 1 II 1 II l 1 1 II 1 1 1 1 II 1 1 1 

Low — 
q =0.0 

; o ; 

O 

^ooc 


i i 1 i i i i 1 i i i i 1 i i i i 1 i ' i i 1 i 

• 

; • • —\ 


1 1 1 1 1 1 1 1 1 1 1 1 1 n 1 1 1 1 1 1 1 1 1 




q Low L :t 

o q & = - 5 

- o * 

o $ ^ ~ 

o 


1 1 | H 1 1 1 1 1 \J i 1 1 1 | 1 1 1 1 1 1 

14ifT>n— T _ 

q„=0.5 

it — a- — ' 

— a 


1 t 1 i H i 1 1 1 1 1 I ! 1 1 1 1 1 t t i i 1 1 


i i 1 h i i i 1 i i i i 1 i i i i 1 i i i i 1 i 



0.1 0.15 0.2 0.25 0.3 0.1 0.15 0.2 0.25 0.3 



log (1 + z) log (1+2) 



Fig. 18. — BCG / magnitude residuals for low-L^ clusters (left panels) and high-L^ clusters (right panels). Residuals are defined as 
Am = mgcQ — m mo d e i, for a passive evolution model with Zf orm = 5 (solid line) assuming f2 m = 0.0, = 0.0 (upper panels) and 
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